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For many-electron systems, we consider a nonequilibrium state (NES) that is driven by a 
pump field(s), which is either an optical field or a longitudinal electric field. For the differential 
optical conductivity describing the differential response of the NES to a probe optical field, 
we derive exact sum rules and asymptotic behaviors, which open wide possibilities for experi- 
ments. In deriving these results, we have also derived universal properties of general differential 
response functions of time-dependent NESs of general systems. 
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Introduction - The optical conductivity tensor cr^^u;) 
describes the response of an equilibrium state to a probe 
optical field. It gives much information on electronic 
properties of condensed matter. 1_12 ) In particular, it 
has been shown that the integrals of Recr^(a;) and 
a/Lm cr^ (uf) over the frequency to are directly related 
to basic properties of the system such as the single- 
particle distribution and band dispersion. 1_12 ) Such rela- 
tions, called sum rules, are therefore useful for exploring 
electron systems, 1-6 ) and have been successfully utilized 
for analyzing a large variety of electron systems. 6 ~ 11 - 1 
However, since an equilibrium state (of each system) is 
uniquely determined by a small number of parameters 
(such as temperature) , the number of controllable param- 
eters that affect the sum (integral) values is very small. 
This fact has severely limited the usage of sum rules. 

This limitation can be removed by considering the op- 
tical conductivity of a nonequilibrium state (NES). A 
NES can be created and driven by a pump field A, which 
is assumed to be an optical field and/or a longitudinal 
electric field (generated by, say, a battery) . The response 
of the NES to a probe optical field is characterized by 
the differential optical conductivity tensor o~ A p [defined 
by eqs. (l)-(3)]. 11-15 ) Unlike equilibrium states and cr^, 
the NES and a A g depend strongly on the magnitude and 
functional form of A(t). Therefore, by tuning A(t) as a 
new controllable parameter, one will be able to make the 
sum rules for o~ A p much more informative than those for 
a ap- However, the problem was that, until now, the sum 
rules for a A g were unknown. 

Note that two different configurations are possible in 
experiments on o A a\ (i) A(t) is turned off before a(t) is 
applied and (ii) A(t) is present when a(t) is applied. We 
here call both configurations pump-probe experiments. In 
configuration (i), the NES (created by A(t) beforehand) 
might sometimes be approximated as a quasi-equilibrium 
state (QES), and the sum rules of are often substi- 
tuted for those of cr^g. 11 -* However, in general, the tran- 
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sient NES is not well approximated as a QES, and this 
substitution has not been justified. In configuration (ii), 
such substitution is obviously wrong because the NES 
driven by A(t) is far from quasi-equilibrium because, for 
example, strong mixing phenomena such as frequency 
mixing take place. Therefore, until now, reliable sum 
rules for o~ A p were unknown in either configuration. 

In this paper, we derive sum rules for a A p(uj) [eqs. (23) 
and (24)], and its asymptotic behaviors [eqs. (25) and 
(26)], for a general class of models for many-electron sys- 
tems. They hold rigorously in both configurations (i) and 
(ii), even when many-body interactions are strong. 

Differential optical conductivity of NESs - Suppose 
that an optical field, described by a vector potential A(t) 
(in the Coulomb gauge), and/or a longitudinal electric 
field, described by a scalar potential </>(r, t), is applied to 
an electron system. Since A and <f> induce optical excita- 
tion and electrical conduction, respectively, the system 
becomes a NES, whose density operator is denoted by 
p A {t). We therefore call A = (A, (f>) the pump field. It 
can be strong such that perturbation expansion in powers 
of A breaks dowi. 13-15 -* Furthermore, we do not assume 
any specific functional form (such as periodicity) for the 
time dependence of A. 

One can study properties of a NES created by A by 
measuring the response to another optical field a(t), 
which we call a probe field. It brings the system into an- 
other NES, p A+a (t). We are interested in the change, 
induced by a(t), in the current density j, 

Aj(t) = {j) A + a -{j) A , (I) 

where (-)f +a = Tr[p A+a (t) ■] and (-)f = Tr[p A (<)-]. 
When a(t) is weak, Aj(t) is well described in terms of 
the differential optical conductivity tensor o~ A p as 

4?«W=£ f t a A p (t-t';t)f (t')dt' + O (f). (2) 
J -°° 

Here, f(t) = —a{t) is the probe electric field, and 
a, ft — x, y, z. Since the NES varies as a function of time, 
so does o A a- That is, a A g depends not only on the time 
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delay r = t — t' but also on t. Furthermore, as eq. (15) 
shows, cr^p is generally a nonlinear functional of A [while 
it is independent of a]. Throughout this paper, the su- 
perscript A, such as those in a\ and (-)f, denotes such 
a functional dependence. Equations (1) and (2) and the 
causality, 



a^(r;t)=0 for r < 0, 



(3) 



define the differential optical conductivity tensor of the 
NES driven by A. It contains much more information 
than that of equilibrium states, o\f^(ui), as we will discuss 
later. 

Experimentally, a(t) is usually taken as monochro- 
matic, and thus f(t) = fe~ lult + c.c. Then, eq. (2) reads 

Aja(t) = JX/jto*)/^-*"* + c.c. + o(f), (4) 

where a^(u>;t) = f°° QQ o A fi (T;t)e lulT dr is the Fourier 
transform (FT) with respect to the time delay r. 16 ' One 
can measure a^g(oj;t) directly by such experiments us- 
ing eq. (27). Since a^p(r;t) is real, Rea^ (uj; t) and 
Im(7^(w;t) are even and odd functions of oj, respec- 
tively. We study sum rules for them. For example, we 
consider 



/oo 
Rea^(w;t)dw, 
-CO 



(5) 



which is called the optical spectral weight. This quan- 
tity is of central interest in many theories and experi- 
ments. 2 ~ 12 ) 

Model and definitions - We consider a many-electron 
system in the presence of electron-electron and electron- 
phonon interactions as well as random potentials. The 
electrons move on a regular lattice, whose dimensionality 
and symmetries are arbitrary. 

We assume that the system is described, in the energy 
scale of interest, by the general Hamiltonian; 



H — H e + H e i + H ee + H ep + Hp . 



Here, H e is the kinetic-energy term of electrons; H e = 
J2 k a e(k)hka, where e(k) denotes the energy disper- 
sion of the band of interest, and hk„ = c ka Ck a - Here, 
Cka = J2i e ik ' 1 cicr / \/N , where cj CT annihilates an electron 
on site I with spin a, and N is the number of unit cells. 
H e i = ^2 1 fj uihia is a random potential (with a random 

on-site energy ui and hi a = c\ ci„), which may be pro- 
duced, for example, by impurities. Furthermore, H ee is 
the sum of electron-electron interactions. We assume that 
H ee is a function of n J(T 's. H ep is the electron-phonon 
interaction, and H p denotes the Hamiltonian of free 
phonons. This general model includes many models such 
as the Hubbard model (for which H ee = U J2 t hi^iii^, 
H e i = H ep = H p = 0). Our results hold irrespective of 
the details and magnitudes of H ee ,H e i and H ep . 

For later use, we define the velocity vector and inverse 
mass tensor as 



h 2 dk a dkj3 



e(k). 



To consider interactions with A and a, we assume that 
the spatial variations of A and a can be neglected. This 
approximation is good in most experimental configura- 
tions. The directions of A, V</> and a are arbitrary. Un- 
der these conditions, we may incorporate the interactions 
with A and a by the Peierls substitution, and the interac- 
tion with <j> by the Coulomb interaction with the charge 
of electrons. Then, the Hamiltonian in the presence of 
A, <j> and a is given by 

fyA+a = Y / e(k- (e/h)A(t) - (e/h)a(t))h ka + H ei 

fe.cr 

+e E ( E - ^ Z < *) + Hee + Hep + Hp. (7) 

I cr 

Here, e is the electron charge, and — en^ s is a background 
charge on site I. By differentiating H A+a with A + a, we 
obtain the current density as 

£ v »( k - ( e / h ) A (t) - (e/h)a{t))h k a (8) 



k,a 



j v a +3a+o(a), 



where 

Jo 



= ^^v a (k - (e/h)A(t))n ka , 



(9) 
(10) 



e 
V 



2 



E m »K k - {e/h)A{t))hkaa {t). (11) 

fe,<7,/3 



When A = 0, j™ represents the diamagnetic current 
induced by a. 2 ~ 6 ) When A^O, the diamagnetic current 
is induced by both A and a, and thus is included in both 
t and j™. 

Since j™ is 0(a), Aj(t) defined by eq. (1) is given by 



Aj(t) = Aj»(t) 



■j m (t)+o(a). 



(12) 



Here, Af a (t) = (f a )f+ a - (j v a )t and j%(t) = 



(6) -J2^d^(t)ap(t), where 



O*) 4E<^ fc -( e / s ) A W)^)t- ( 13 ) 



fe.cr 



For a simple cubic lattice, for example, Ea^a^W 1S 
proportional to the expectation value of the kinetic en- 
ergy. 

While j m (t) responds to a(t) instantaneously, Aj v (t) 
responds with a finite delay as 

A f a (t) =E fs%{t-t';t)ap(t')dt' + o(a). (14) 

J -°° 

Here, $^(T;i) is the response function describing the 
differential response of Aj v (t) to a(t). We denote its 
FT with respect to the time delay r by E^g(w; t). Since 
f(t) = -a(t), eqs. (2) and (12)-(14) yield the differential 
optical conductivity tensor as 



v a p(u;t) 



u + iO 



Z*Ju;t)-d*M) 



(15) 



Both E^g and d^g are nonlinear functionals of A, and so 
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Universal properties of response functions of time- 
dependent NESs - To derive sum rules for a A p, we note 
that E^g in eq. (15) should satisfy all the universal prop- 
erties that were found in ref. 13 for general response func- 
tions of general systems. Since ref. 13 assumed steady 
NESs driven by a static pump field, we here generalize its 
theory to time-dependent NESs, which are realized, for 
example, by the application of a time-dependent pump 
field. For this general discussion, we omit vector and ten- 
sor indices. 

We denote the pump and probe fields by A(t) and 
a(t), respectively. In nonequilibrium statistical mechan- 
ics (e.g., in the Kubo formula 1 ) and in refs. 13-15), it is 
usually assumed (implicitly) that an observable of inter- 
est is independent of a(t). However, we here consider the 
general case where an observable of interest, denoted by 
Q a it), is a function of a(t), because this is the case for j a 
given by eq. (8). Then, by expanding Q a (t) m powers of 
a(t), we obtain 

Qa(t) =Q + Qia(t) + o(a), (16) 

where Q and Q\ are operators independent of a{t). We 
have obtained such an expansion in eq. (9), where Q = 
f a and Q\a(t) = j™. The response to a(t), AQ a ^ = 

(Qa(t))f +a - (Qa(t))t, is therefore given by 

AQ a{t) = AQ(t) + (Q^fait) + o(a), (17) 

where AQ(t) = (Q)f +a - (Q)f- Since the response func- 
tion of the second term on the right-hand side is sim- 
ply given by {Qi)f 1 let us consider the non-trivial term 
AQ(t). Unlike (Qi)f, AQ(t) depends on p A+a {t) (the 
NES in the presence of both A and a). We therefore have 
to use the theory of ref. 13 to evaluate AQ(t). 

When a(t) is sufficiently weak, AQ(t) responds to a(t) 
linearly as 

AQ(t)= I ^ A (t-t';t)a(t')dt' + o(a). (18) 

J — oc 

This and the causality condition, <& a (t; t) = for r < 0, 
define the differential response function $ A (r;t) of the 
NES. Its FT with respect to the time delay r is denoted 
by E a (uj; t). It is straightforward to generalize the theory 
of ref. 13 to the case where A and the NES are time- 
dependent. We then obtain the following results. 
The dispersion relations, such as 

-4- Im» A ( W ';t) — , (19) 

-oo w w 77 

are satisfied. Furthermore, the sum rules 

/OO 7 
ReE A (u;t)— = (C)t, (20) 
- oo 1" 

j QcIm~>;t) - {C) A } *± = (D)f (21) 

hold. Here, C = [R,Q]/ih and D = —[Q, [R, H A + 
H']]/H 2 , where R denotes the operator that couples to 
a(t) via the interaction term —Ra(t), H A is the Hamil- 



tonian of the target system in the presence of A [such 
as eq. (7) with a = 0], and H 1 is the interaction be- 
tween the target system and other systems such as heat 
reservoirs and electric leads. 13 ) In general, these opera- 
tors (such as Q and R) are additive operators or their 
densities. 13 ' 15 ^ Equation (21) also gives the asymptotic 
behavior for large oj as 

LjImZ A (uj;t) -> (C) A . (22) 

In deriving these results following ref. 13, we have used 
the von Neumann equation for the density operator of a 
huge system, which includes not only the target system of 
interest but also environments and a source of the pump 
field, as well as all interactions among them. [Although 
such a huge system is analyzed, we have successfully de- 
rived, as in ref. 13, the relations among quantities of only 
the target system.] Therefore, these results are rigorous 
and apply to all physical systems, as long as the linear 
relation given by eq. (18) holds. 13 ' 14 ) 

Main results - Let us apply the above results to o~ A p 
of the system described by eq. (7). By expanding H A+a 
in powers of a(t), we find that R = Vjp for ap{t). For 
E A fj , which is the FT of $^ of eq. (14), Q = j*. The 
sum rules for a A (i are obtained from the properties of 
Z A through eq. (15). 

For the optical spectral weight [defined by eq. (5)], 
eq. (19) for uj = yields 17 ) 

W A f} {t)=*d%{t). (23) 

Note that this result relies only on eqs. (15) and (19). 
That is, this sum rule is derived only from the causality 
[eq. (3)] and the specific form of the current [eq. (9)]: No 
other relations are necessary for deriving this sum rule. 
For co Im a A p, on the other hand, eqs. (15) and (20) yield 
the following sum rule: 

/OO 
{wIm^(w;t)-d^(t)}d W = 0. (24) 
-oo 

This and eq. (22), respectively, give the asymptotic be- 
haviors for large to as 

wlm<7^(w;t) d%{t), (25) 

uj 2 Rea A fj {uj;t) -> 0. (26) 

Equations (23)-(26) are our main results. They are 
rigorous (to the same degree as the Kubo formula is) 
within the general model denned by eq. (7), even when 
A(t),(j)(t),H ee ,H e p and H ei are strong. For example, our 
results hold for any possible phases of the system that is 
described by eq. (7). That is, our results are completely 
valid as long as the target system is well described by 
the Hamiltonian of eq. (7). Conversely, if experimental 
results disagree with our results, it means that the sys- 
tem is not described by eq. (7) (because, say, transition 
to another band takes place). Such rigor seems impor- 
tant for the application of the sum rules and asymptotic 
behaviors. 

Note that the effects of H ee , H ep , H ei and <f) on the sum 
and asymptotic values appear only through the distribu- 
tion function {hk a ) A . In contrast, the effects of A on 
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the sum and asymptotic values appear not only through 
{nka)t but also through m~^(k—(e/h)A). In either case, 
the decoherence of electrons affects the sum and asymp- 
totic values only through the broadening of (hka)f- 
Possible applications - For (7°^, the sum rule for 

VfQ ee /^Re^M^ reads VQ/ir = <Q ee 

(e 2 /^)E fc , CT ^(fc)(^)° q - 1_6) For each system, <Q 
depends only on the temperature T and doping density 
n<j. In pump-probe experiments, in contrast, W^/n = 
dap(t) can De studied as a function of T, n<j and A. This 
opens wide possibilities for studying many-electron sys- 
tems. For example, suppose that an ordered phase is re- 
alized as an equilibrium state. By measuring cr°^, one 
obtains the value of for the ordered phase. Then, a 
static A — (0, 4>) is applied to induce a DC electric cur- 
rent while keeping T equal to that for A = (by, for 
example, using a good heat sink). By measuring afi^ 
by applying a(t), one now obtains, from eq. (23) or (25), 
the value of $f = (e 2 /V) J2 k a m^(fe)(n fc(T )( » for 
a non-ordered phase, because the order would be de- 
stroyed by the electric current if \V(j>\ was larger than 
a certain value. One thus obtains the values of d a p with 
and without the order at the same T and rid- Alterna- 
tively, suppose that no order is present in an equilibrium 
state. Then, a coherent optical field A = (A(t), 0) is ap- 
plied. This would induce an electron-hole (eh) correla- 
tion. Hence, by measuring cr£g, one obtains the value of 
d£p for the state with the eh correlation. 

Method of measuring a^g(co;t) - a£g(w;t) can be 
measured, for example, by the following process. 

Step 1: Prepare the system in some initial state at 
an initial time t = 0. Apply a pump field A(t) only, 
and measure the current density j(t) continuously for a 
sufficiently long time. Then, turn off A(t), and at another 
initial time prepare the system in the same initial state as 
that at t = 0. Redefine the origin of time (t = 0) as this 
new initial time. Apply the same pump field A(t) again, 
and measure the current density j(t) continuously. By 
repeating these procedures sufficiently many times, one 
obtains many independent records of j(t). The average 
of these records gives 

Step 2: Perform the same sequence of experiments us- 
ing the pump and probe fields instead of the pump field. 
Here, the pump field A(t) is taken to be the same as that 
of Step 1. One then obtains (j)f +a . From this and the 
result of Step 1, one obtains Aj(t) = (j)f +a - 
If one takes the probe field cLS cL monochromatic one, 
f(t) = fe~ lut + ex., and if one takes / parallel to the 
/3-axis (i.e., f a = fS a0 ), then eq. (4) yields Aj a (t) = 

Step 3: Perform the same sequence of experiments us- 
ing the same pump field and another (phase shifted) 
probe field f'(t) = f'e~ iut + ex., where /' = if. One 
then obtains (j)f +a . From this and the result of Step 
1, one obtains Aj'(i) ee (j)f +a ' - According to 

cq. (4), it is expressed as Aj' a (t) = a^ g (nj;t)ife~ lult + 
cx. + o(f). 

From these experimental results, one can evaluate 



^ (<*>;*) usin S 

a aB (uj, t) = hm 2j— t . (27) 

Concluding remarks - Our results hold in both 
configurations (i) and (ii), which were discussed in 
the introduction. In configuration (i), cq. (23) reads 
W&(t) = (^e 2 /V)j: k rT m^(k)(n ka )t. Comparing this 
with the corresponding result for cr^g, 1-6 ) W^g = 
(7re 2 /y) J^ fe cr m~^(fe)(n fc(T ) cq , we find that the result 
for W\{t) is obtained simply by replacing the equi- 
librium electron distribution (nfc cr ) eq with the nonequi- 
librium one (fika)t- Hence, the analysis of the pump- 
probe experiments in ref. 11, which substituted the sum 
rule of W c a l for that of W^(t), is now justified. In 
configuration (ii), on the other hand, cq. (23) reads 
W*(t) = (7re 2 /V)£ fe;(7 m-^fc - (e/h)A(t))(h k „)f. 
Since the pump field enters the inverse mass tensor, the 
simple replacement of (rifc (T ) cq with (fik<j}f in the sum 
rule of W^i does not yield the correct result. 

Finally, we point out that the present results can be 
generalized. Suppose that the current density takes a 
general form; 

J a = J*-Y J t) %a (} {t)+ {a). 

Here, an d are arbitrary vector and tensor oper- 
ators, respectively, which may be functions of A. [Equa- 
tion (9) takes this form.] Then the sum rules eqs. (23) 
and (24) are respectively generalized as 

/oo 
Rea^(cj;t)du J = 7T(D^)t, (28) 
-OO 

rLlma^(u;;t)-(D^)t}doj = 0. (29) 

J — oo ^ J 

Furthermore, generalizations to the case where the probe 
field is a longitudinal AC electric field and to higher-order 
responses [following ref. 15] are straightforward. 
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